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Abstract 

Two families of solutions of a generalized non-Abelian Toda lattice are considered. These 
^ solutions are expressed in terms of quasideterminants, constructed by means of Darboux and 

binary Darboux transformations. As an example of the application of these solutions, we con- 
I ^ I sider the 2-periodic reduction to a matrix sine-Gordon equation. In particular, we investigate 

the interaction properties of polarized kink solutions. 

> 

ON 1 Introduction 
in 

There has been great interest in noncommutative versions of some well-known soliton equations, 
such as the KP equation, the KdV equation and the Hirota-Miwa equation [1-12]. Often, these 
noncommutative versions are obtained simply by removing the assumption that the coefficients in 
the Lax pair of the commutative equation commute. 
J> The non-Abelian Toda lattice 

^ f/n,x + ^nVn+X " VnVn = 0, (1.1) 

Vn.t + t/n-1 - = 0, (1.2) 

was ffi'st studied in [13]. A Darboux transformations for this system was given by [14]. In [15], the 
following generalization 

C/n,x + f/„K+l - VnUn = 0, (1.3) 

K,t + anf/n-1 - f/nOn+l = 0, (1.4) 

was studied and the Darboux and binary Darboux transformations were obtained. We note that in 
general, an is not a scalar, but is independent of t. In the case that J7n, I4i and an are scalars, it is 



easy to show, by setting Q„+if7„ = e " and eliminating V^, that (1.3|-(1.4) becomes the standard 
two-dimensional Toda lattice equation 



'n,xt 



2e"^" - e"^"+i = 0. (1.5) 
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Introducing new variables Xn where 

Un = ^nX^^i, Vn = X^^xX^^ , (1-6) 



(1.3|-(1.4| can be rewritten as 

[Xn^xX^ + 0.nXn~lX^ ^ — XnX^^-^^ttn+l = 0. (1.7) 

From now on, we will refer to (|1.7|) as the non-Abelian Toda lattice. One type of quasideterminant 



solutions of (1.7) were found in [16]. We will show how these (quasiwronskian) solutions arise 
from the Darboux transformation and consider a second type of quasideterminant, which we call 
quasigrammian, solutions obtained using the binary Darboux transformation. 

It is well known that the 2-periodic reduction of the standard two dimensional Toda lattice leads 
to the scalar sine-Gordon equation. In the same way, the 2-periodic reduction of the non-Abelian 



Toda lattice (1.7) leads to a noncommutative sine-Gordon equation. This equation has been studied 
already in a number of papers [16-23] concerning both the matrix and the Moyal product versions. 
Here we only consider in detail the matrix version. 

Recently, a matrix KdV equation was considered in [24] . A multisoliton solution was found by 
using the inverse scattering method. In particular, the properties of one- and two-soliton solutions 
expressed in terms of projection matrices were investigated. We will apply some of these ideas to 
the matrix sine-Gordon equation to study the interaction of its kink solutions. 

The paper is organized as follows. In Section 2, some properties of quasideterminants used in 
the paper are described. In Section 3, we present quasiwronskian solutions to the non-Abelian Toda 
lattice constructed by iterating Darboux transformations and in Section 4, we present quasigram- 
mian solutions to the system by using the related binary Darboux transformation. In the rest of 
the paper we consider the 2-periodic reduction to a matrix sine-Gordon equation. In particular we 
consider the matrix kink solutions obtained from the quasigrammian solutions, we show that kink 
solutions for the matrix sine-Gordon equation emerge intact from interaction apart from change of 
polarization and phase. 



2 Preliminaries 

In this short section we recall some of the key elementary properties of quasideterminants. The 
reader is referred to the original papers [16,25,26] for a more detailed and general treatment. 



2.1 Quasideterminants 

An n X n matrix A over a ring TZ (noncommutative, in general) has quasideterminants written 
as \ A\ij for i,j = 1, . . . ,n, which are also elements of TZ. They are defined recursively by 



1^1 



A' 



(2.1) 



In the above rj represents the ith row of A with the jth element removed, c*- the jth column with 
the ith element removed and A^'^ the submatrix obtained by removing the ith row and the jth 
column from A. Quasideterminants can be also denoted as shown below by boxing the entry about 
which the expansion is made 



1^1 



^,3 
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The case n = 1 is rather trivial; let A = (a), say, and then there is one quasideterminant 

a 

c d , 



1^1 1,1 = = a- For n = 2, let A = ( ^ ) , then there are four quasideterminants 



l^li.i 

\A\2,l 
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b 





d 



a — bd ^c, |^|i,2 
c - db~^a, \A\2,2 ■ 
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d 



b — ac ^d, 
d — ca~^b. 



Note that if the entries in A commmute, the above becomes the familiar formula for the inverse of 
a 2 X 2 matrix with entries expressed as ratios of determinants. Indeed this is true for any size of 
square matrix; if the entries in A commute then 



1^1 



^,3 



. w+j det(^) 
^ ' det(^^J) ■ 



(2.2) 



In this paper we will consider only quasideterminants that are expanded about a term in the 
last column, most usually the last entry. For a block matrix 

A B 
C d 

where d £ TZ, A is a square matrix over TZ of arbitrary size and B, C are column and row vectors 
over IZ of compatible lengths, we have 



A B 



C d 



d-CA-^B. 



2.2 Noncommutative Jacobi Identity 

There is a quasideterminant version of Jacobi's identity for determinants, called the noncommuta- 
tive Sylvester's Theorem by Gelfand and Retakh [25] . The simplest version of this identity is given 

by 

^ B C An A n A u A n 

(2.3) 
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we have the homological relation 
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(2.4) 



2.3 Quasi-Pliicker coordinates 

Given an {n + k) xn matrix A, denote the ith row of A by Ai, the submatrix of A having rows with 
indices in a subset / of {1, 2, ... , n + k} by Aj and ^{i,...,n+fe}\{i} by A^. Given i,j £ {1,2, . . . ,n + k} 
and / such that = n — 1 and j ^ /, one defines the (right) quasi-Pliicker coordinates 



rUA) :-- 



Ai 




Ai 


Ai 


ns 


A, 



-1 



Ai 

Ai W 



A, 1 



(2.5) 



for any column index s G {1, . . . , n}. The final equality in (2.5 ) comes from an identity of the form 



(2.3) and proves that the definition is independent of the choice of s. 
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Remark. A useful consequence of (2.5) is the identity 



A' 





-1 


A' 





A' 







A' 


1 


A^ 


1 




A^ 






(2.6) 



which shows that quasideterminants of this form may be inverted very simply. 



3 Solutions obtained by Darboux transformations 

The non-Abelian Toda lattice (|1.3|)-(|1.4|) has Lax pair 



t>n,x = VnCpn + a„(/>n-l, 
i>n,t = Un(l)n+l- 



(3.1) 
(3.2) 



Let On,i ,i = 1, ...,N be a particular set of eigenfunctions of the linear system and introduce the 
notation 0„ = {9n,i,-'' ^Gn,N)- The Darboux transformation, determined by particular solution 



9n, for the non-Abelian Toda lattice is 



^1 



Vn = Vn + anOn-lOfi^ — OnOn+l'^n+l, 

This may be iterated by defining 

(t)n[k + 1] = Mk] - en[k]en+l[k]~^(l)n+l[k]. 
Xn[k + 1] = en[k]en+l[k]-^Xn+l[k], 



where (/>„[!] = (/>„,X„,[1] = Xn and 



In particular, 



4>n^] — 4>n — 9n,l0n+l,l^n+l, 
Xn[2] = On,lOn+l,l^n+l- 



(3.3) 
(3.4) 
(3.5) 
(3.6) 



(3.7) 
(3.8) 



(3.9) 

(3.10) 
(3.11) 



In what follows, we will show by induction that the results of N repeated Darboux transforma- 
tion (f)n[N + 1] and Xn[N + 1] can be expressed as in closed form as quasideterminants 



e 



n+N (Pn+N 



Xn[N + l] = (-1 



e„ [0 
e„+i 



e 



n+N 



X. 



n+N- 



(3.12) 
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The initial case = 1 follows directly from ( 



Jacobi identity (2.3) and the homological relation (2.4) we have 



3.10 )-( 3.11 ). Also using the noncommutative 



bn[N + 2] = 0„[iV + 1] - 9n[N + l]en+l[N + 1]" Vn+l[iV + 1] 

©n+l 



en 

©n+l (pn+l 

&n+N 4>n+N 

e 

©n+l (Pn+l 



©n. Sn,N+l 
©n+l ^n+l,Af+l 



©n 



+2 



7n+l,N+l 
%i+2,N+l 



n+N ^n+N,N+l ©n.+AT+l On+N+l,N+l 



© 



n+l 



(/'n+l 



©n+2 (j^n+i 



©n+AT+l '/'n+AT+l 



'n 'rn 



©n ^n,Af+l 
©n+l ^'n+l,Af+l 

©n+Af ()n+N,N+l 



©n+Af '/'n+AT 

©n ^n,Af+l 
©n+l ^n+l,Af+l 4>n+l 



©n+Af+1 ^n+Af+l,Af+l 4>n+N+l 



©n+l 
©n+2 



7rt+l,iV+l 
'^n+2,N+l 



©n+AT+l &n+N+l,N+l 



©n+l 
©n+2 

71+Af + l 



(Pn+l 
4>n+2 



© 



»n+Af+l 



and 



Xn[N + 2] = 0„[7V + l]en+i[N + l\-^Xn+i[N + 1] 

©n 



(-1) 



N 



7n,N+l 



©n+l ^n+l,Af+l 



+1 



7n+l,Ar+l 



(-1) 



N 



7n,N+l 



©n+l ^n+l,Af+l 



©n+AT ^n+Af,Ar+l 



©rt+1 

©n.+2 ^n+2,Af+l 

n+Af+1 ^n+Ar+l,Af+l 

©n.+l ^n+l.AT+l 



© 



©n+l 



©n+2 
©n+Af+1 1 



n+A^+l 



©n+AT ^n+Af,Ar+l 
©n+7V+l ^n+Af+l,Ar+l 



-1 



X, 



n+Af+1 



and then using the quasi-Pliicker coordinate formula (2.5), we get 



(-1) 



Af+l 



©n ^'n,Ar+l L0_ 

©n+l ^n+l,Af+l 



X, 



n+Af+1- 



©n+AT+l ^n+Af+l,Af+l 1 

This proves the inductive step and the proof is complete. 



4 Solutions obtained by binary Darboux transformation 



The linear equations (3.1) and ( |3.2[ ) have the formal adjoints 

-1pn,x = V^lpn + aj^+iVn+l, 



n-l- 



(4.1) 
(4.2) 
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Following the standard construction of a binary Darboux transformation, one introduces a 
potential n„ = fi((?!>„, ipn) satisfying the three conditions 



^{<f>n,'^n)t = -1pliUn(f>n+l, 
fin - = -tpli(l>n- 

A binary Darboux transformation is then defined by 

+ = MN] - en[N]n{en[N],pn[N])-'n{^n[N],pn[N]), 
^N + i] = v„[7V] - pn[N]n{en-i[N],pn-i[N])-''n{en-i[N],^n-i[N])\ 

Xn[N+l] = {I + en[N]n{en[N],Pn[N])-^Pn[N]'')Xn[N], 

where 0„[1] = ^„,V'n[l] = V^n, Xn[l] = and 



(4.3) 
(4.4) 
(4.5) 



(4.6) 
(4.7) 
(4.8) 



(4.9) 



Using the notation G„ = (^n,i, • ■ ■ , ^n,iv) and P„ = {pn,!, ■ ■ ■ , Pn,N), it is easy to prove by induction 
that for iV > 1, 



Q(e„,p„) n(0„,p„) 



j)(e„_i,p„_i)t j^(e„_i,Vn-i)^ 



and 



We may thus after N binary Darboux transformations we obtain 

^(©nj Pn) Pil 



Xn 



(4.10) 
(4.11) 

(4.12) 
(4.13) 



In fact, we can prove the above results by induction. 

Xn[N + 2] = (/ + en[N + l]n{en[N + 1], Pn[A^ + l])-'Pn[iV + l]^)Xn[N + 1] 



^(©n^Pn) ^{Gn,N+l,Pn 



7n,N+l 



^(©ri) Pfi 



f^(^n,Af+l, Pn 



^(©n-l, Pn-l) 
^^(0n-l,/>n-l,Ar+l) 



P' 



yn,N+l 



^{&n, Pn,N+l) ^{On,N+l, Pn,N+l) 
^(®n) Pn) Pn 



On 
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Noticing 



^{&n-l, Pn-l) 
^{&n-l, Pn-l,N+l) 



n 

Pn,N+l 



~iPn,N+l 



n{en-i,Pn-i,N+iMen-i,Pn~i)~'Pl){i + enn{en,Pn)-'Pl) 

pIi,N+1 + Pn-l,N+l)^{Qn-l, Pn-l)'^Pl 

{fl{Qn, Pn,N+l )-n{en-l,Pn-l,N+l)Men,Pn)-'Pl 

j)(G„_i,p„_i,iv+i)^^(e„_i,p„-i)-i(i^(e„_i,p„_i)-n(e„,p„))n(e„,p„)-^pt 

plN+l+^{en,Pn,N+lMQn,Pn)-^Pl 

^{en,Pn) 
^{&n, Pn,N+l) 



Pn,N+l 



we have 



Xn[N + 2] 



5 Pn) Pn 
^{Qn, Pn,N+l) ^{On,N+l, Pn,N+l) Pn^N+1 
©n &n,N+l —I 



5 Matrix sine-Gordon equation and its kink solutions 

It is well known in the commutative case that one may obtain reductions by imposing periodic 



conditions on the On- Similarly in non-Abelian case, one can make periodic reductions of (1.7). 



From now on, we only consider the case that X„ is a d x d matrix and a„ = Idxd and so (1.7) is 

{Xn,xXn + Xn-lXj^ ^ — XnX^^-^ = 0. 



(5.1) 



The simplest such reduction has period 2, that is, we take X„ 
system 

(Xq^xXq + XiXq ^ — XqX-^ ^ = 0, 
{Xi^xX^^)t + XqX^-^ — XiXq^ = 0. 



Xn+2 and (5.1) gives the 



(5.2) 
(5.3) 



We call this a non-Abelian sinh-Gordon equation since in the commutative case, it will be seen that 
Xq = X^^ = Fi/Fq and then 9 = 21og(Fi/i<o) satisfies the standard sinh-Gordon equation 



^xt 



4 sinh 6. 

By changing 9 i9, we can also obtain the sine-Gordon equation 

9xt = 4 sin 61. 



In what follows, we will construct solutions to (5.2)-(5.3) by reduction of the solutions (4.13 



of the non-Abelian Toda lattice (5.1). It is clear that (5.1) has vacuum solution Xn = I and (4.13 
gives the quasigrammian solutions 



Xn 



^(©rii Pn 



(5.4) 
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where 9n,i and satisfy 



(5.5) 



and is defined by (4.3|-(4.5). We choose the simplest non-trivial solutions of (5.5) 



Pn,i — -^iPi 6 



where Ai and Bj are d x d matrices and then we obtain 



A^-BjPi [ p. 



Qj-Pi \Qj 



The choice of constant of integration as 5ijl is needed to effect the periodic reduction we will 
shortly make. This can also be written as 

^ (I)" (...(|)%f!fe'"--^*")- 

Now using the invariance of a quasideterminant to scaling of its rows and columns (see e.g. [26]), 
we get 



Xr, 



I H — ^-^e 

Qj-Pi 



'gi\2 — ... — / QN \2 
'Pi' ^Pn' 

Pi = —Qi = Xi for i = 1, • • • ,N. Therefore, the non-Abelian sinh-Gordon equation has the solutions 



It is obvious from this expression for Xn that it is 2 periodic when i^)'^ = ■ ■ ■ = (^)^ = 1, that is. 



Xn 



X.x-j-UT 



{Bje ^/) 



(5.6) 



(5.7) 



Prom now on we will assume that Ai = I are real and Bj = irjPj, where rj are real scalars, are 
pure imaginary matrices. In this case, it follows that Xq and Xi are complex conjugate to one 
another. For this reason we introduce 



X — Xn — X^ 



J-' (A.+Aj)'' 

{B,e'^'"~^. 



(5.8) 



Next, we will derive matrix kink solutions for the matrix sine-Gordon equation using the method 
applied to study the soliton solutions of the matrix KdV equation in [24] . To get a visual repre- 
sentation of the solution we will consider the matrix W{x, t) defined by 



iW^ = X^X- 



XxX 



(5.9) 



We choose this dependent variable so that in the scalar case W = 0, the solution of the sine-Gordon 
equation. 
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For TV = 1, gives 



X = I + B { I - ^e- 



(5.10) 



We first assume further tliat P is a projection matrix (i.e. satisfies = P)- Tliis clioice allows us 
to calculate the inverse matrices in the above expression explicitly using the formula 



(I-aP)-i =1 + 

where a 7^ 1 is a scalar and P is any projection matrix. 
In this way we find that 



aP 

1^' 



(5.11) 



X = I + 



irP 



and hence 



^2Xx+ft 



4XP 



ir/2 



cosh(2A(x + t/A2 -0)) 



where (j) = log(r/2)/2A. Note also that XX = I. 

Taking one final step, we integrate to obtain the one-kink solution to the matrix sine-Gordon 
equation 



W = 4Parctan(e 



2A(z+t/A2- 



(5.12) 



Remark. For the one-kink solution (5.12), we call the projection matrix P its polarization and (p its 
phase. In the scalar case, if we choose P = 1, (5.12) is simply the one-kink solution to the standard 
sine-Gordon equation. 



For N = 2, expanding X by the definition (2.1 ), we can rewrite X as 



X = I + {Bie ^1 ,B2e 



"^^'^ (A. + A,)' 



1 / -Xi'x-^t- 



2x2 \ e 



2-1 



1 1 / -Xix-j-t J 

/ + [Lie ^1 , L2e ^2 ) 



-X2X—^t 



>-2 I 



= 1 + L1+L2, 

and hence 



Li(e ^1 /- -Si) 

J I 2X2X^^1 J 1 „ s 



Ai -h A2 
Ai 



L2-B1 = B\, 
L1B2 = B2. 



Ai + A2 

In the expressions Bj = irjPj, j = 1,2 we assume that Pj are the rank-1 projection matrices 



Pi 



Pj ® qj _ Pjij 
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and the d-vectors pj and qj satisfy the condition {pj, qj) ^ 0, we can solve for Li and L2 by using 
(5.11 ) to obtain 

Li = ^^i^^(A2i?2 + (Ai + X2)92l)Bu 

a 

L2 = i^i±M(Aii?i + (Ai + X2)giI)B2, 
9 



where 



and 



where 



Therefore 



5 = (Ai + A2) gig2 + Xi\2rir2a, where a 



{Pi,Qi){P2,q2) 



5,=e2V^-^ for J = 1,2, 



1 

X + — 



X = I+ i^i±M(AiSiS2 + A2S2S1 + (Ai + X2){giB2 + g2Bi)). (5.13) 

g 

We now investigate the behaviour of X as t ^ ±00. We will use the fact that W is invariant 
under the transformation X XC for any constant matrix C and assume, without loss of gen- 
erality, that < Ai < A2. In the calculations that follow, we will demonstrate that kinks emerge 
from the interaction and undergo phase-shifts as in the scalar case. In addition however, we will 
see that there are changes of polarization, in other words, amplitudes may also change as a result 
of the interaction. 

First we fix 9i. Then 62 = O1 + (I/A2 — 1/Af)t and so as t — > —00, 



gi e2Aiei _ ^ 

As t ^ -|-oo, using the invariance of W , we obtain 

2B2i (Ai + A2)(2AiBiS2 + 2X2B2B1 - (Ai + A2)r2Sii) - 4aXiX2riB2i 



Xr.I + 



r2 r2{Xi + X2Ygii - 2aXiX2rir2 



^ _ (Ai + A2)(2Ai^i^2 + 2X2B2B1 - (Ai + A2)r2^ii) - 4aAiA2ri^2i ^ 2i^\ \ (j^ 252Z 



where 



r2(Ai -I- X2Ygii - 2aXiX2rir2 \ r2 J J \ r2 

iriPi 



^ 2 



nipi,qi) £, pi^qi 

2X2 {pi,q2) . 2A2 {P2,qi) 

pi=pi-T\ — ^^'2' qi = qi-j\ — r^? v^2- 

(Ai + A2) (^2,92) (Al + A2) (P2,g2) 
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This shows that 

W ~ 4Pi arctan(e^^i(^i-'^^)), t -oo, 
W ~ 4Pi arctan(e2^i(^i-'^^)), i ^ +00 

where 

1 , n ,+ 1 , ri{pi,qi) 

Similarly, fixing 62, we have 



X I -\ — t ^ —00, 



2 



where 



Then 



Bo iroPo 



(P2,92) (P2,g2) 

2Ai (P2,gi) . 2Ai (pi,g2) 

(Ai + A2) (pi, gi) (Ai + A2) (pi, gi) 



~ 4P2 arctan(e2^2(e2-</.2 )) ^ t ^ -00 
W ~ 4P2 arctan(e2^2(e2-<^f))^ ^ ^ ^^^^ 



where 



2X2 2{p2,q2) 2A2 2 

The above calculations show that W{x,t) decomposes into the sum of two kink solutions as 

.j. The phase shifts Aj = (f>'^ — (f>J 



t ±00 and the jth kink solution propagates with the speed 1/A^. The phase shifts Aj = (f>^ — • 



for the kink solutions are 
where 

4AiA2a 
^" (Ai + A2)2- 

Remark. In a similar way to the matrix KdV equation in [24], we find that the matrix amplitude 
of the first kink solution changes from 4Pi to 4Pi and the matrix amplitude of the other one 

changes from 4P2 to 4P2 as t changes from -co to +00. If (^1,^2) = (P2P1 = 0) or {p2,(li) = 
(P1P2 = 0), there is no phase-shift, however the amplitudes may change. In the case that both 
P1P2 = P2P1 = 0, there is neither phase-shift nor change in amplitude and so the kink solutions 
have trivial interaction. 

To illustrate the above, we will consider the case d = 2, i.e. the 2x2 matrix sine-Gordon 
equation. We choose Ai = 1, A2 = 2, ri = r2 = 1, and 



Pi 



1 "M p ^ 1 2 -2 

r ^ 3 V -1 1 
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Pi 




Figure 1: Asymptotic forms for kink 1 





-1 1 



P.' 



Figure 2: Asymptotic forms for kink 2 




The analysis above shows that P^ = Pi 



= P2 and 



-1 1 
-4 4 




-1 1 



For convenience, rather than plotting the kink W given by (5.9), we plot the derivative Wx and 



refer to it as a soliton. In Figure [T] the asymptotic forms of the matrix soliton 1 are plotted as 
t — > ±00. The first plot exhibits the amplitudes given by P^ and the second, those of P^ . Similarly, 
in Figure [2] we show the same plots for soliton 2. 



6 Conclusions 

In this paper, we have considered a generalized non-Abelian Toda lattice and presented quasiwron- 
skian and quasigrammian solutions obtained means of by Darboux transformations and binary Dar- 
boux transformations respectively. Then we imposed a 2-periodic reduction on the non-Abelian 
Toda lattice to derive a noncommutative sine-Gordon equation. By using a method similar to 
that developed in [24] for the matrix KdV equation, we obtained kink solutions for the matrix 
sine-Gordon equation from the quasigrammian solutions of the non-Abelian Toda lattice. Then we 
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investigated the interaction properties of matrix kink solutions. It is known [27] that the change 

of matrix amplitude of solitons for the matrix KdV equation gives rise to a Yang-Baxter map. It 
would be interesting to investigate whether there is a similar result for the matrix sine-Gordon 
equation. 
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